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INTRODUCTION 
In  the  present  stage  of  development  of  dynamic  meteor- 
ology, the  study  of  the  large-scale  nearly-horizontal  raotionB 
in  the  atmosphere  appears  most  susceptible  to  analytical  (and 
computational)  treatments   In  this  paper  we  avoid  the  study 
of  larger  scale  phenomena  which  entails  Icnowledge  of  the 
physical  ( thermodynamic al)  nature  of  the  general  circulation 
of  the  atmosphere  and  the  study  of  smaller  scale  meteorological 
motions  involving  viscosity  and  turbulence.   Now,  two  of  the 
most  well-established  observations  in  meteorology  are  1)  the 
atmosphere  is  in  a  state  of  nearly-hydrostatic  equilibriimi 
and  2)  predominantly  its  motion  is  nearly-geos trophic,  except 
in  a  relatively  thin  viscous  layer  on  the  earth's  surface. 
The  first  observation  means  that  the  pressure  at  most  points 
in  the  atmosphere  arises  substantially  from  the  vxeight  of  the 
vertical  colixran  (unit  area)  of  air  above  the  point  and  that 
the  acceleration  term  in  the  vertical  (with  respect  to  the 
earth's  surface)  component  of  the  momentum  equation  can  be 
neglected  to  the  lovjest  order  of  approximation.   This  is 
related  to  the  observations  that  the  effective  depth 
(wherein  lies  most  of  the  gaseous  mass)  of  the  atmosphere 
is  small  compared  to  the  earth's  radius  and  that  its  motion 
is  predominantly  horizontal.   The  second  observation  means 
that  the  horizontal  pressure  gradient  is  nearly  balanced  by 
the  horizontal  component  of  the  apparent  Coriolis  force^due 
to  the  earth's  rotation  and  that  the  acceleration  terms  in 
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the  horizontal  components  of  the  momentum  equation  can  be 
neglected  to  the  lowest  order  of  approximation.,  VJe  shall 
see  that  this  is  consistent  x^rith  the  additional  physical  ob- 
servation that  the  motion  is  slow  compared  to  any  other  char- 
acteristic velocity,  for  example,  the  soiond  or  gravity  wave 
velocity;  that  is,  the  only  characteristic  velocity  of  any 
consequence  is  the  particle  flow  velocity.  By  these  approx- 
imations, the  possible  time-dependent  motions  are  reduced 
considerably 0 

Recently  several  meteorologists,  aided  by  some  ideas 
and  suggestions  of  Sossby  (19^^-0)  have  made  important  con- 
tributions in  attempting  to  obtain  sim.pler  formulations  of 
atmospheric  motion  by  introducing  these  (and  other)  physical 
approximations  into  the  complicated  equations  of  gas  dynamics, 
Charney  (19l!.8),  Thompson  (1952)  and  others  have  discussed 
rather  thoroughly  the  orders  of  magnitude  of  the  physical 
variables  and  have  obtained  approximate  formulations  v;hich 
are  being  subjected  to  the  high-speed  numerical  methods 
required  for  short-range  weather  forecasting  in  the  middle 
latitudes.   One  of  our  prirnary  raotives  in  this  paper  is  to 
demonstrate  that  approximate  formulations  can  be  obtained  by 
a  consistent,  systematic  procedure.   This  procedure  auto- 
matically yields  higher  order  approximations  without  requiring 
additional  physical  assumptions <, 

In  part  I  of  this  paper  we  derive  approximate  flovr  equa- 
tions incorporating  the  hydrostatic,  geostrophic  and  low- 
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velocity  approximations  vrith  the  simplest  possible  model  which 
we  believe  to  be  a  physically-meaningful  description  of  some 
atmospheric  (and  oceanic)  motions.   The  method  which  i-ie   use 
is  apparently  a  rather  general  and  powerful  procedure  of 
systematically  introducing  physical  approximations  into  a 
well-formulated  but  complicated  mathematical  description  of 
a  physical  system.   The  approximate  systems  which  result  are 
asymptotic  representations  in  some  sense  depending  on  the 
particular  problem  considered.   The  method  consists  of 
combining  two  formal  mathematical  procedures:   1)  transforming 
each  variable  (dependent  and  independent)  by  multiplying  it 
by  g.   sm^all  pararaeter  raised  to  an  arbitrary  exponent  and 
2)  making  a  power  series  expansion  xjith  respect  to  that 
parameter  of  each  of  the  dependent  variables.   By  introducing 
a  sufficient  n'-omber  of  consistent  physical  approximations 
the  exponents  can  be  evaluated.   The  second  procedure,  by 
itself,  is  the  widely-used  linear  perturbation  method  of 
approximation  which  always  retains  the  underlying  local 
geometrical  properties  of  the  original  sjrstem..   The  first 
procedure,  by  stretching  or  contracting  each  variable,  allovjs 
an  ordering  of  the  terms,  dictated  by  the  physical  approxima- 
tions; this  is,  essentially,  Prandtl's  boundary  layer 
procedure,  x^rhiGh  implies  that  we  are  seeking  certain  similarity- 
type  solutions.   Hoviever,  the  small  parameter  need  not  occur 
explicitly  in  the  original  differential  equation  system.   The 
combined  ordering  approximation  procedure  is  particularly 


suitable  for  those  physical  systems  which  are  described  by 
complicated  differential  equations  but  in  which  the  orders 
of  magnitude  of  the  terms  for  the  particular  phenomenon 
of  interest  are  kno;m.   This  is  a  weak,  but  often  natural 
and  sometimes  sufficient,  condition  on  the  properties  of 
the  desired  approximate  solution  so  that  the  differential 
equation  system  which  it  satisfies  can  be  derived.   Sub- 
stantially, this  method  (if  not  completely  the  viewpoint) 
is  the  same  as  that  used  by  Friedrichs  (19ii-3)  and  Keller 
(I9I4.8)  in  their  study  of  surface  wave  phenomena.   Another 
example  of  the  ordering  approximation  procedure  is  given 
by  Morikawa  (195^4-)  •   The  present  paper  is  its  sequel  in 
method  and  viewpoint;  in  both  of  these  physical  problems, 
the  retention  of  a  non-linear  formulation  to  the  lov;est 
order  of  approximation  is  apparently  essential  for  express- 
ing fully  the  implications  of  the  physical  approximations. 
A  related  approximation  procedure,  in  which  both  the  depen- 
dent and  independent  variables  are  expanded  with  respect  to 
a  small  parameter,  has  recently  been  exploited  successfully 
by  Lighthill  (19^-^9)  and  others. 

In  Part  II  of  this  paper  we  study  the  physical  system 
described  by  the  loi/est  order  of  approximation.   This  system 
is  governed  by  a  physically-simple  conservation  lav;  of 
geostrophic  motion,  although  the  non-linear  mathematical 
formulation  presents  difficulties.   Hovjever,  this  conserva- 
tion law  suggests  a  precise  concept  of  atmospheric  vortices 


5. 

in  a  rotational  flow,  in  close  analogy  with  ordinary  point 
vortices  in  a  tvjo- dimensional,  incompressible,  irrotational 
tlo\-J,      Consequently,  we  are  able  to  describe  the  motion  of 
these  atmospheric  point  vortices  by  a  system  of  first  order 
ordinary  non-linear  differential  equations.   Visualized  in 
terms  of  atmospheric  vortices,  the  large-scale  nearly- 
horizontal  motion  of  closed  high  and  low  pressure  systems 
may  be  conceptually  simplified.   Recently,  James  (1950) 
has  proposed  a  different  concept  of  atmospheric  vortices  on 
a  more  phj/^sical  basis. 


PART  I 

1,  Formulation, 

In  a  completely  systematic  derivation  of  approximate 
meteorological  equations,  one  would  begin  with  the  mathe- 
matical formulation  of  atmospheric  motion  on  a  spherical 
earth  based  on  the  equations  of  gas  dynaraicso"^  But  this 
would  involve  additional  physical  approximations,  primarily 
thermodynamical  ones,  which  are  not  so  xifell-established  as 
the  hydrostatic  and  geostrophic  approximations.   So  v;e  begin 
by  making  overall  physical  assumptions  which  are  believed  to 
be  valid  in  establishing  a  primitive  approximate  model.  These 
are  1)  making  a  tangent  plane  approximation  to  flow  on  a 
spherical  earth  and   2)  considering  a  mechanical  model  by 
assijming  that  thermodynamical  effects  are  secondary.   These 
assianptions  have  precedents  in  the  work  of  L.  Prandtl,. 

We  consider,  then,  an  inviscid,  homogeneous  atmosphere 
(with  a  free  surface)  on  a  horizontal  plane.   Friedrichs  (I9I4-8) 
has  shown  how  the  hydrostatic  approximation  can  be  Introduced 
into  such  a  physical  ^stem  (but  non-rotating),  yielding  the 
equat  ions  of  the  well-known  shallov;  water  theory.  Therefore 
we  relegate  to  Appendix  I-A  a  more  obvious,  simpler,  deriva- 
tion.  Introducing  the  horizontal  component  of  the  Coriolis 
force,  the  resulting  equations  of  nearly-horizontal  free- 
surface  flow  represent  the  simplest  possible  physical  model 
of  large-scale  horizontal  motions  in  the  atmosphere.   This 
physical  system"^ is  described  by  the  equations  of  continuity 
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atid  momentuin  foi?  the  depth  h(x,3*-^t)   and  the  horizontal 

velocity  components  u(x>y>t)  and  v(x,y>t)  in  the  X  and 
y  directions^  respectively^ 

f+h(u^+Vy)=0  (la) 

fH  +  gh^  .  fv  =  0  (2a) 

^  +  gh^  +  f  u  =  0  (2b) 


where 


dt 


^=  (  )t  -^-u  .  (  ),  -H  V  M  )^,  (2c) 


means  differentiation  along  particle  paths  in   (x,y,t)  space, 
g  is  the  acceleration  of  gravity  and  f  =  20.sin  <j>  is  the 
Coriolis  parameter  in  terms  of  the  earth's  angular  velocity 
n  ahd  the  latitude  ahgle  (j)  measured  from  the  equators   f 
is  fixed  at  the  desired  point  of  tangency  to  the  earth's 
surface i   Subscripts  denote  partial  dif f erentiat ion*   In 
addition  the  pressure  p  and  vertical  velocity  w  vary  linearly 
in  the  vertical  direction  z* 

p  =  gy)(h  -  z)  (3a) 

where  p  is  the  density  (constant)  and 

w  =  -z(u^  +  V  )4  (3b) 

Equations  (la)  and  (2a, b)  are  a  system  of  three  first  order 
non-linear  partial  differential  equations  for   (h,u,v)   of 
hyperbolic  type.  That  is,  the  characteristic  manifold^ 
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j^(x,y,t)  =  constant,  along  which  certain  discontinuities  are 
propagated,  is  defined  by 

The  first  term  represents  the  particle  paths  in   (x,y,t)  space 
and  the  quantity  within  the  bracket  describes  the  Monge  cone 

(cf.  Courant  and  Hilbert).   The  characteristic  gravity  wave 

1/2 
velocity   (gh)  '    is  of  the  same  order  of  magnitude  as  the 

velocity  of  sound  in  the  lower  atmosphere.   Upon  this  property, 
partly,  depends  the  appropriateness  of  the  equations  (la)  and 
(2a, b)  in  describing  atmospheric  motions.   Usually,  the  con- 
tinuity equation  (la)  is  replaced  by  the  vorticity  equation 
obtained  by  differentiating  equation  (2a)  with  respect  to  y 
and  equation  (2b)  with  respect  to  x,  taking  the  difference 
between  them,  and  eliminating  the  divergence   (u  +  v  ) 

^   y 

between  the  resulting  equation  and  equation  (la).   Thus  we 
obtain 

ft  (H^)  =  °  (^^) 

where  ^  =  v  -  u   is  the  vertical  component  of  vorticity. 
The  system  of  equations  (lb)  and  (2a, b)  is  more  directly 
descriptive  of  atmospheric  motions  than,  but  not  essentially 
different  from,  that  of  equations  (la)  and  (2a,b),   Since 
either  system  presents  formidable  analytical  and  computa- 
tional difficulties,  we  are  forced  to  consider  additional, 
reasonable  physical  approximations. 
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2.   The  Geostrophlc  and  Loxj-Velocity  Approximations, 

We  now  introduce  the  geostrophic  and  low-velocity  approx- 
imations into  the  system  of  equations  (la  or  lb)  and  (2a,b) 
by  the  ordering  approximation  procedure  (cf.  Introduction), 
In  terms  of  the  stretching  and  expansion  parameter  a,  the 
geostrophic  approximation  states  that  the  first  term  (actually, 
three  terms),  in  each  of  equations  (2a)  and  (2b),  is  of  higher 
order  in  a   (loi^/er  order  of  magnitude)  than  the  second  term 
which  is  of  the  same  order  of  magnitude  as  the  third  term  (to 
the  lowest  order  of  approximation).   The  low-velocity  approx- 
imation states  that  we  retain  the  particle  derivative 
d(  )/dt  intact  and  make  a  perturbation  expansion  on  the  atmos- 
phere at  rest,  i^e,,  on  the  solution 

(h,u,v)  =  (h^,0,0).  (5) 

We  note  that  this  is  an  exact  solution  of  equation  (la)  and 
(2a, b)  and  also  of  equations  (lb)  and  (2a, b)  if  f  =  constant. 
The  evaluation  of  the  stretching  parameter  exponents  for 
each  variable  is  carried  out  in  Appendix  I-B,   The  simple  ■ 
result  is  that  the  time   t   is  the  only  stretched  variable. 
This  means  that  the  combiiBd  geostrophic  and  low-velocity 
approximations  imply  an  asymptotic  representation  for  large 
time.   That  is,  to  the  lovrest  order  of  approximation,  the 
flow,  which  would  initially  move  in  the  direction  of  negative 
pressure  gradient,  is  directed  normal  to  this  direction  by 
the  horizontal  Coriolis  force  when  the  description  is  given 
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in  terms  of  the  appropriate  time  scale  t.   Thus,  the  time   is 
transforraed  into 

t;  =  at  (6a) 

and  the  system  of  equations  (la  or  lb)  and  (2a, b)  is  trans- 
formed by  replacing  d(  )/dt  by  d(  )/dT:  v;here 

^  =  a(  )^  +  u(  )^  +  v(  )y.  (6b) 

The  solution   (h,u,v)   and  the  characteristic  solution  |_  are 
expanded: 

h(x,y,a-)  =  h^  +  ah^^^(x,y,T:)  +  a2h^^^x,y,'r)  +  ...  (7a) 

u(x,y,T)  =  au^^^(x,y,T)  +  a^u^^^ (x,y,T)  +  ...  (7b) 

v(x,y,T)  =  av^-'-^x,y,T)  +  a^v^^^  (x,y,T)  +  ...  (7c) 

|(x,y,'c)  =Io  ■"  af^^^x^y.-r)  +  a2$^^^x,y,'c)  +  ...  (7d) 

We  carry  out  the  approximation  by  putting  equations  (7a,b,c) 
into  the  transformed  system  of  equations  (la  or  lb)  and  (2a,b) 
and  eq;jating  coefficients  of  llke-povxers  of  a  in  the  usual 
way.   We  do  the  same  with  equation  (7d)  and  the  transformed 
characteristic  equation  corresponding  to  equation  (1|)« 

The  solution  (h^-^  ^  jU^"'- ^ ,  v^""-^ ) ,  the  lowest  order  (first 
order  in  a)  approximation,  satisfies 
(1)    f■^^        .u(l) 


d 


(^(1)  -  £i2_ —  )  =  0,   f  =  constant^,   (8a) 


dT   ''^        h 

o 


gh^l^  -  fv^^^  =  0  (8b) 

gh^^^  +  fu^^^  =  0  (8c) 
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where 

dl^.  (  )^.u<l'  .  (  )^-v(l'  .  (  )y         (8d) 

and  the  vorticity  ^^-^^   =  v^'   -  ^i^    •     Since  f  =  constant, 
equation  (Sb,c)  gives 

K^^^   =  i  A  h^^^  (8e) 

where  the  Laplacian  A  =  (   ^xx  "^  ^  ^ty'   ^^^  *^®  divergence 

u^^^  +  v^^^  =  0.  (Sf) 

X     y 

^      ".sfies 


The  characteristic  solution  (t)^    satl 


dT 

Equation  (8a)  is  the  non-linear  conservation  equation  con- 
sistent xifith  the  geostrophic  equations  (8b,c),   The  character- 
istic equation  (8g)  states  that  only  the  particle  paths  remain 
as  true  characteristics  since  the  Monge  cone  has  degenerated 
into  a  plane  t  =  constant;  this  is  the  precise  statement  of 
the  low-velocity  approximation.   The  resti'iction  that 
f  =  constant   is  apparently  a  consequence  of  the  tangent  plane 
approximation  (cf.  Section  1).   The  only  consistent  way  which 
we  can  see  to  relax  this  condition  is  to  consider  a  better 
approximation  to  flow  over  a  spherical  earth. 

The  solution   (h^'^  ,  u^   ,  v^    ' ) ,    the  second  order  (  in  a) 
approximation,  satisfies 
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iil  „U)  .  4!!,  ,  ,„U)|.  ,  .(2)|.H,'l'  -  4^)  (9a, 

O  /-■  X  O 

^(1)  2  .(1)  h  ?'^' 

o  h 


(2)  (2)        d^-'-^u^-'-^ 

(1)„(1) 


gh^'^^  -  fv^^^  +S_^  =  o  (9b) 


g,(2),,,(2)  ,^llL^  =  o.  (90) 

(?)     (2)  (2) 

Prom  equation  (9b, c)  the  vorticlty  ^^  ''  =  v^   "  v    becomes 

^(2)  =£^h(2)  -  2J(u(l^v(l^)       (9d) 
where  the  Jacoblan  J(a,b)  =  a  b  -  %^x  ^^*^  ^^^®  divergence 

„(2) .  ,(2)  = .  i  liiisiii . .  I.  liivii  ^ .  !4;^.(,^, 

^x      y      f    dT       h     dT        h  '^^^^ 
•^  o  o 

The  vertical  velocity  to  this  order  of  approximation  is 

w'2'  =-z(u(2'  .v<2')  =^'.  (ff) 

^  X         X  h  vv   / 

o 

Equations  (fa,b,c)  give  a  system  of  non-homogeneous  linear 
differential  equations  with  non-constant  coefficients  in 
terras  of  the  lowest  order  solution   (h   ,  u^  ,    v    )•  VJe 
do  not  consider  this  non-p:eostrophic  system  further  in  this 
paper  although  it  may  be  of  some  interest  to  meteorologists^^ 
vrho  wish  to  consider  higher  order  approximations. 
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Appendix  I-A,  Derivation  of  the  Non-Linear  Shallow  Water 
Equations,^ 

Consider  the  system  of  equations  describing  three- 
dimensional,  incompressible  free-surface  flow  over  a  hori- 
zontal plane,  z  =  Og   The  continuity  and  momentum  equations 
satisfied  by  the  velocity  components   (u,v,w)   and  pressure 
p  =  pP   (constant  density  p)  are 

u__  +  v_  +  w_  =  0  (I-l) 

X    y    z 

-±-   +  P_  =  0  (I-2a) 

dt     X 

d^v   _ 

-±-   +  P_  =  0  (I-2b) 

dt     y 

d-w   _ 

-J-   +  p_  +  g  =  0  {I-2c) 

dt     z 

where 

d.(  ) 

-\^      =  (   )_  +  u  •  (   )_  +  V  •  (   }_  +  w  .  (   )_,  .  (I-2d) 
t  X  y  z 

The  boundary  conditions  at  the  free   ("z  =  h{x,y,t ) )   and  fixed 
(z  =  0)   surfaces  are 


P(x,y,h,t)  =  0  (I-3a) 

w(x,y,o,t)  =  0  (I-3b) 


w(x,y,h,t)  =  ^  (I-3c) 

dt 


where 


••■.■•Ki 


Ik. 

^^^   =  (   )_  +  u(x,7,h,t)  •  (   )_  +  v(x,y,h,t)  •  (   )_.  (I-3d) 
dt       t  X  y 

In  addition  we  take  the  horizontal  components  of  vorticity  to 
be  zero: 

u_  -  w_  =  0  il-l\.a) 

Z      X 

v_  -  w_  =  0,  (I-ij-b) 

z    y 

The  shalloi^r  water  equations  are  derived  by  a  stretching  trans- 
formation only ; V 

X  =  e^x,  y  =  e  y,   z  =  e^z,      t  =  e  t  (1-5) 

u  =  e^u,  V  =  e%,  w  =  e%,  P  =  e^^P 

vjhere  the  linbarred  quantities  are  the  stretched  variables, 
8  is  a  small  parameter,  and  the  exponents  are  constants  to 
be  deten:nined  by  the  physical  approximations.   Transforming 
equation  (I-l),  (1-2),  (1-3)  and  (l-k)    into  the  stretched 
variables  by  equation  (1-5)*  the  physical  approximations  which 
we  make  are : 
1)   Particle  paths  are  preserved,   3y  equation  (I-2d) 

dt  ^  ^^  7 

+  e^-'^  •  w  .  (  ) 

z 

yields 

d  =  a  +^  =b  +m  =  c  +n.  (I-6a) 
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2)  Hydrostatic  pressure  approximation  •   By  equation  (I-2c) 

e^-  .  tf  .  ,0.,  ,^  ,  g  .  „ 

yields 

c+q=0,      d+n>0,  (I-6b) 

3)  The  x*-  and  y-xnomentum  equations  are  preserved.   By  equations 
(I-2a)  and  (I-2b) 


and 


dt    "    "x 


d+m   ^31  J.  b+o  T,  _  ^ 
^    •  dt  ■"  ^   "  ^y  -  ° 


yield 


d  +/  =  a+q,   d+m  =  b+q,  (I-6c) 


Equations  (I-6a,b,c)  are  six  equations  for  eight  unknowns. 
For  example,  we  can  solve  in  terms  of  n  and  q: 

a  =  b  =  n  -  3q/2,   c  =  -q,   d  =  n  -  q 

(I-7a) 
/'  =  m  =  q/2. 

The  inequality  equation  (I-6b)  gives 

2n  -  q  >  0.  (I-7b) 

Thus  the  simplest  stretching  transformation  satisfying  equation 
(I-7a,b)  is  given  by 
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q  =  0,   n  =  1.  (I-7c) 

Then  the  shallow  water  equations  result  from  transforming 
equation  (I-l),  (1-2),  (1-3)  and  (I-/|)  by 

X  =  ex,  y  =  ey,   z  =  z,  t  =  et 

(1-8) 
u  =  u,   V  =  V,   w  =  ew,   P  =  P 

and  setting  e  =  0  in  the  transformed  equations.   The  irrota- 
tlonallty  conditions  (I-I|.a,b)  specify  that  the  velocity  com- 
ponents  u  and   V  are  independent  of   z  to  this  order  of 
approximation.   Summarizing 

^  +  h(u  +  V  )  =  0  (I-9a) 

dt     X    y 

f  +gh^  =0    ,  (I.9b) 

g4-gh^  =  0  (I-9C) 

P  =  g(h  -  z)  (I-9d) 

w  =  -z(u^  +  V.  ).  {I-9e) 

Equations  (l-9a,b,c)  are  the  well-knox^m  shallow  water  equations 
where   (h,u,v)   are  functions  of   (x,y,t)   only  (cf,  Friedrichs 
19k-Q   and  Keller  19ii.8),  VJe  note  that  the  physical  interpre- 
tation of  the  parameter  e  is  still  free. 
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Appendix  I-B.  The  Geostrophic-Loiir  Velocity  Stretching 
Transformation, 

For  the  meteorological  equations  (lb)  and  (2a, b)  consider 
the  stretching  transformation   (h,u,v)  (original  )v  —c» 
(h,u,v)  (stretched)   and   (x,y,t )  (original )  -— >  (^jVIjT)  (stretched), 
with  respect  to  a  small  parameter  a 


(I-IO) 


,.    a        b         c  ■ 

£,  =  a  X,     K|  =  ay,   T  =  at 

—    /    —    m    r-    n, 
u  =  au,   v  =  av,   h  =  an 

combined  with  the  perturbation  expansion  of  (h,u,v) 

h  =  h^  +  ah^^^(C,r|,'i:)  +  a^h^2^4,y|,T:)  +  ... 

u  =      au'^^^,n,T)  +  a^u^^^C,)n,T)  +  ...         (I-ll) 

V  =      av^^U^,y^,'^)    +   a2v^^^(^,r|,T)  +  ...   , 

Then  we  demand  that   (h^-'-,u^-'-^,v^-'"^)   satisfy  the  following 
conditions : 

1)  Particle  paths  are  preserved. 

a^  .  1^  =  a°  •  (  ),  .  a-^-^^^  •  u^^^  .  (  )^ 

+  ^b.m.l  .,(!).() 

1 

yields 

c=a+>^  +  l=b+m  +  l.  (I-12a) 


»  ' 
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2)  Vortlclty  Is  preserved. 

yields 

a  +  m  =  b  +  /  .  (I-12b) 

3 )  Geostrophic  approxlmatlon» 


and 


d-u 
yields 

a+n=in,  c+/>in  (I-12c ) 

and 

b   +  n  =   /,  c   +  m  >    /  ,  (I-12d) 

Equations    (I-12a,b,c ,d)    are   four  independent   equations   for 
six  unlcnowns;    solving  in  terms   of     c      and     n 

a  =  b   =  ^(c   -  n  -   1),     i  =  m  =  ^{c   +  n  -   1).  (I-13a) 

The  inequalities,  equations  (I-12c,d),  reduce  to 

c  >  0,  (I-13b) 

The  simplest  stretching  transformation  satisfying  equations 
(I-13a,b)  is  given  by 
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n  =  0,  c  =  1  (I-13c) 

which  makes   a=b=^  =  in  =  0,  ¥e  obtain  the  extremely  simple 
result  that  time  is  the  only  variable  stretched,  the  trans- 
formation being  t  =  at. 

We  note  that  if  v/e  attempt  to  incorporate  the  above 
approximations  1),  2)  and  3)  in  the  meteorological  equations 
by  stretching  transformation  only,  that  is,  in  the  same  way 
that  we  obtained  the  shallow  water  equations  in  Appendix  I-A, 
we  end  up  with  two  alternatives:  i)  either  f  =  constant  and 
h  =  0  or  ii)  h  progresses  without  change  of  shape  in  the 
east-west  direction  with  a  velocity  gh(f~  )   which  is  an 
order  of  magnitude  larger  than  the  gravity  wave  velocity 

1/2 

(gh)  '  ,   ii)  is  physically  unacceptable,  but  1}   is  consistent 

with  our  results,  i.Co,  the  atmosphere  at  rest,  (h  ,0,0),   For 
completeness,  the  stretching  transformation  for  this  case  is 
defined  by  the  following  relationships  bet^^:een  exponents: 

a  =  b  =  -^(c  -  n),   ^  =  m  =  ^(c  +  n); 

and  the  inequality: 

c  >  0,  (I-li^b) 

For  example  we  can  choose 

n  =  0,   c  =  1  (I-lii-c) 

;^rhich  makes   a=b=/=m  =  l/2   (cf,  eqns.  (I-13a,b,c  )  )  , 
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PART  II 

3,   The  Geostrophic  Conservation  Equation, 

Although  the  immediate  analytical  difficulties  of  a  non- 
linear approximate  formulation  are  evident,  an  apparent  ad- 
vantage is  the  concise  physical  statement,  equation  (8a),  which 
results.   Consequently,  a  mathematical  advantage,  which  may  not 
be  apparent  a  priori,  is  that  the  basic  structure  of  the  desired 
solution  is  essentially  retained,  ^^rhile  the  undesired  pro- 
perties are  pushed  into  the  backgroiind. 

The  conservation  equation  (8a)  can  be  regarded  as  a  single 
third  order  non-linear  differential  equation  for  the  perturbed 
depth  h^   ,  eliminating  u   ,  v^     and   J^     \jj   equations 
(8b,c,d,e);  or  introducing 

11/  =  I  h^^^  (10) 

the  system  of  equations  (8a,b,c)  becomes 

(A  -  ^  )^  =  0,  (11a) 


dT 


u^^^  =  "Vl/y  (lib) 

v^^^  =  vl/^  (lie) 

where   k^  =  f^gh^.   We  call  i/\   -  k^)   the  Helmholtzian 

operator  and  equation  (8a  or  11a),  the  geostrophic  conservation 

equation.   The  geostrophic  equations  (llb,c)  show  that  ^ix,j,x) 

AD 
is  the  stream  function   ( — ^^   =  v|/  )  with  the  identity  given  by 

equation  (10),   In  particular,  we  note  that  the  sign  (plus 
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or  minus)  of   v]/  and  the  perturbed  depth  h^     is  the  same. 
In  addition,  we  point  out  that  the  flox-J  described  by  equations 
(8)  or  (11)  is  divergence-freeX/  (see  equation  (8f))  which 
implies  that  the  vertical  velocity  w^    -  0   (by  equation  3b), 
Hence,  there  is  no  interchange  of  potential  and  kinetic  energy 
to  this  order  of  approximation;  this  interchange  is  described 
by  the  second  order  approximate  equations  (9)» 

Since,  in  the  middle  latitudes,  there  visually  exist  pre- 
vailing westerly  winds,  it  may  be  convenient  (but  not  necessary) 
in  using  numerical  metliods  to  replace   \I/  by  [\!/  (y)  +  ^jf] 
where   \1;  (y)   is  prescribed  by  the  loioxvn  westerly  flox^r 
U(y)  =  -(\!i  )  •   Then  equation  (11a)  becomes 


||-.,„.„(i),|...a)|_j(^.K2„ 


(lid) 


-  (U   -  K.^U)\1/  =  0, 

yy   ^  '^x 


From  the  viewpoint  of  studying  the  validity  of  the 
geostrophic  conservation  equation  in  describing  atmospheric 
motions,  we  discuss  three  possible  approaches:   1)  numerical 
method,  2)  linearization,  and  3)  atmospheric  point  vortices. 

1)  Numerical  Methodc^'^  A  problem  of  primary  interest  to 
meteorologists  is  that  of  vjeather  prediction.   For  example 
from  weather  data  at  a  given  time  over  a  region,  say  the 
United  States,  we  wish  to  predict  its  motion  there  for  some 
later  time.   If  i-re  desire  a  prediction  over  the  entire  region 
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we  must  also  give  some  data  continuously  (t  >  0)   at  the  bound- 
aries of  the  region.   Thus,  we  consider  initial  and  initial- 
boundary  value  type  problems  for  the  equations  (lla,b,c)  by 
finite  differences.   The  conservation  form  of  equation  (11a) 
suggests  what  are  reasonable  conditions  (data)  and  a  compu- 
tational procedure.   For  the  initial  value  problem  in  the  entire 
x-y  plane  we  prescribe  ^{x,j,0)      (or  h^   (x,y,0))o  Then  for 
T  =  0,  \-ie   can  calculate  at  each  grid  point  the  velocity 
(u^  f^        ),  giving  us  the  initial  slope  of  the  particle  path 
in   (Xjy^T)   space,  and  {/\  -   k  )^,  which  is  conserved  along 
the  particle  path  according  to  equation  (11a) e  We  can  now 
calculate  the  distribution  P(x,y),  say,  of  {/\   -  K  )\]/  at 
grid  points  after  a  time  step  /\^  %,    limited  by  /^t  <  /\x/u 

and  /\T  <   AyA^ax'  ^^®^^   (A  x,Ay)   is  the  grid  spacing 
and   (^^^^>v^^)   is  given  by  the  initial  distribution  of 

(u^-'-^,v^-'-M.   If  v!/(x,y,A'^)   is  prescribed  at   (x,y)  =  co 
we  can  determine  it  in  the  entire   (x,y)   plane  for  t  =  A'^ 

as  the  solution  of   (A  -    )^  =  F(x,y),  and  so  on  for  the 
next  time  step.  For  a  finite  region  (initial-boundary  value 
problem),  we  must  prescribe  both   vl/(x,y,T:)   and  the  vorticity 
A^(x,y,T).  on  the  boundar^^  for  all  t  >  0,  where  A^  is 
prescribed  only  at  those  points  at  wMch  the  velocity 
(u   ,v^   )   is  directed  into  the  region.   The  computational 
procedure  outlined  here  apparently  differs  from  that  presently 


(1) 
max 
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used  by  meteorologists X^ (see  Charney  19^9  and  Thompson  1953 )• 

2)  Linearization.   Some  insight  into  the  nature  of  solutions 
of  the  Geostrophic  Conservation  Equation  (11a)  can  be  obtained 
by  linearization.   For  example,  we  perturb  (h^   >u   \v   )   on 
a  known  west-east  (westerly)  flow  of  velocity  U(y): 

h^-^^(x,y,T)  =  h^Cy)  +  ah^(x,y,T)  +  a^h2(x,y,T)  +  ,..   (12a) 

(1)  p 

u    (x,y,T)  =  U(y)  +  au^(x,y,T)  +  a  U2(x,y,T;)  +  ...    (I2b) 

(1 )  2  '' 

V   '(x,y,i:)  =        av^(x,y,T)  +  a  V2(x,y,T)  +  ,.,  ,  (12c) 

Putting  equations  (12a, b,c)  into  equations  (8a,b,c,d,e ),  the 
first  order  approximate  solution  h, (x,y,T)   satisfies  the 
linear  third  order  partial  differential  equation  with  non- 
constant  coefficients 

(|?  -^  ^  fe^^A  -  ^^^)hi  -  (h-^)   (^  -  K^)U  =  0      (13a) 

with  the  geostrophic  conditions 

U  =  -g(h^)y/f,  u^  =  -g(h^)y/f,  v^  =  g(h^)yf.     (13b) 
For  U  =  constant,  equation  (13a)  reduces  to 

or  written  in  a  slightly  different  form 

ih   -^  u  ^)A\   -  ^^i\)^   =  0  (13d) 

This  is  the  same  differential  equation  (with  f  =  constant) 
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studied  by  Charney  (191-1-9),  Thompson  (1952)  and  others.  However, 
there  are  three  differences:   1)  the  time  scale  is  different, 

2)  the  interpretation  o£    K        is  changed,  and  3)  in  their 
linearization  the  dependent  variable  corresponds  to  a  first 
order  perturbation  whereas  the   h,   in  equation  (13d)  is  a 
second  order  (in  a)  approximation  of  the  depth  h,  cf,    _ 
equations  (7a)  and  (12a)»,  Another  linearization  worth  men- 
tioning is  given  in  Appendix  II-A;  this  is  obtained  by  the 
ordering  approximation  procedure  with  a  perturbation  expansion 
of   (h,u,v)   on  a  given  west-east  flow  starting  directly  with 
the  system  of  equations  (1)  and  (2)  such  that  the  geostrophic 
approximation  (for  the  velocity  component  u,  only)  holds  up 
to  the  first  order  in  a.  A  comparison  of  these  lineariza- 
tions would  be  instructive  but  we  do  not  carry  out  the  analysis 
in  this  paper — partly  because  we  note  that  a  linear  perturba- 
tion on  U  implies  that  the  desired  approximate  flo\-j   is  of 
higher  order  of  magnitude o   Based  on  our  present  knowledge, 

we  believe  that  it  is  preferable  to  regard  the  westerly  flow 
U  to  be  of  the  same  order  of  magnitude  as  the  desired  lowest 
order  approximation,  eeg*,  as  described  by  equation  (lid), 

3)  Atmospheric  Point  Vortices,   The  form  of  the  Geostrophic 
Conservation  Equation  (11a)  suggests  another  approach  by 
analogy  to  two-dimensional  incompressible  rotational  flox-x. 
This  analogous  flo^^^  is  described  by  the  continuity  equation, 

u  +  V  =  0  (li+a) 

X     y 


^  u        ...  ^,       v 
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which  implies  the  existence  of  a  stream  function   \1;  such  that 

the  velocity  components   u  =  -"ij/^  and   v  =  "ij/  :  and  the  vorticity 

J  ^ 

X  =  V     -  u     =  /\lj     satisfies 

ft  An?  =  0  Hkh) 

where  the  particle  derivative 

^=  (  )t  -u  .  (  )^-^v(  )y.  (11,0) 

Equation  (llj-b)  states  that  the  vorticity  is  conserved  along 
particle  paths  in  two-dimensional  incompressible  flow.   The 
analytical  difficulties  of  studying  the  Geostrophic  Conservation 
Equation  (11a)  are  similar  to  those  of  studying  equation  (ll+b). 
However,  in  the  study  of  a  large  class  of  fluid  flow  phenomena, 
experience  has  shown  that  a  reasonable  approximation  is  that 
the  flow  is  irrotational  (^  =  /\^  =   0)  almost  everywhere  except 
in  small  isolated  regions  which  can  often  be  approximated  by 
point  singularities  in  two  dimensions  or  line  singularities  in 
three  dimensions.   Then  this  approximation  t'];ives  a  good 
representation  of  the  physical  picture  except  in  the  immediate 
vicinity  of  the  singularities,  IJell-lmown  examples  are 
Th,  von  Karman's  representation  of  the  "vortex  street"  flow 
behind  a  circular  cylinder  and  L,  Prandtl's  "horseshoe  vortex" 
approximation  of  a  finite-span  airfoil.  We  hypothesize  that 
there  are  atmospheric  flows  which  may  be  regarded  as  almost 


.(1)  _  fh^^Vh^)  =  (A  -  "^ 
=  0  almost  everywhere  except  in  small  isolated  regions 


"Helmholtzian-free";  that  is,  {^^^'    -   fh^-^Vh^)  =  (A  '  '^  )^ 
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representable  by  point  singularities.   The  consequences  of  this 
assvmiption  are  considered  in  the  next  section.   At  this  point  we 
note  that  the  G-eostrophic  Conservation  Equation  (11a)  is  in- 
herently two-dimensional;  that  is,  if  we  try  to  specialize  to 
one  space  dimension  (x  or  y),  equation  (11a)  becomes  time- 
independent. 

Some  remarks  are  appropriate  here  about  the  possible  choice 
of  the  parameter  a  for  the  three  approaches  outlined  above. 
Thus  far  the  only  evident  restriction  on   a  is  that  as 
a  — ^  0,  the  solution   (h,u,v)  — *»  (h  ,0,0),  the  state  of  the 
atmosphere  at  rest.   So  we  have  some  freedom  in  the  physical 
interpretation  of  a  depending  both  on  our  approach  and  on 
the  particular  physical  problem  which  we  wish  to  consider. 
Since   a  is  related  not  only  to  the  physical  parameters 
(g  and  f )  which  appear  as  coefficients  of  the  original  system 
of  differential  equations  (la  or  lb)  and  (2a,b),  the  physical 
interpretation  must  arise  from  the  initial  and/or  boundary 
conditions.   For  1)  numerical  methods  or  2)  linearization,  we 
can  relate   a  to  the  maximum  initial  velocity,  e«g«, 

a  =  ||[u(l^x,y,0)]2  +  [  v(^^x,y,0)]  ^j^^^  (l5a) 

(grad  h^   )max. 
or  another  alternative, 
a  =  (gh^)-l/2  .|[u(l^x,y,0)]2  +  [  v^^  ^x,y,  0)  ]  ^j^^^  .   (i^b) 


I.Vjf 
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For  atmospheric  vortices  (see  equation  (l6a)  and  sequel)  of 
strength  y^  >  i  "=  1»  • ' » >^  ^®  can  take 

(l^c) 


f 
a  =  — 


(y- ) 
'i  max 


S^o 


Since  our  approximation  procedure  is  based  on  the  introduction 
of  a  characteristic  vertical  depth  h   (and  not  a  character- 
istic horizontal  length),  equation  (l5b)  is  a  more  natural 
ciioice  than  equation  (l5a)»   This  preference  can  be  illustrated 
by  the  following  discussion:   For  the  purpose  of  argument,  sup- 
pose we  introduce  a  horizontal  length  S  as  well  as  a  charac- 
teristic vertical  length  H   (or  h  )  and  a  characteristic 

o  O 

horizontal  velocity  V,  and  non-dimensionalize  the  momentum 
equation  (2),   Then  the  low-velocity  and  geostrophic  approx- 
imations state  that 

V^  «  gH  =  fSV 

which  can  be  rewritten 

in  <<  1   ^^d   fs  <<  1- 

Since  we  have  introduced  the  depth  h   in  a  rather  natural 
way  (cf.  Section  2,  equations  (5)  and  (7a)),  we  focus  our 
attention  on  the  first  inequality  (small  Froude  number)  Instead 
of  the  second  (small  Rossby  number )e 

At  the  present  stage  of  development  of  the  ideas  in  this 
paper,  vje  prefer  to  regard  both  the  small  parameter   a  and 


'  <  •   ■'■  "  f: 
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the  depth  h   as  physical  parameters  which  are  at  our  disposal 
in  order  to  fit  as  closely  as  we  can  the  mathematical  model  of 
a  homogeneous  atmosphere  with  a  free  surface  to  the  actual 
physical  model  of  nearly-horizontal  atmospheric  motions.   The 
primary  role  of  a  is  to  change  the  tiine  scale  (equation  {6a) ); 
the  primary  role  of  h  ,  from  the  viewpoint  of  the  lowest  order 
approximation  (the  Geostrophic  Conservation  Equation  (8aJ  or 
(11a)),  is  to  relate  the  balance  of  potential  energy  to  kinetic 
energy  (cf,  remarks  follov/ing  equation  (llc))o 


„i .  'f  -t  n  ■ 
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U.,     The  Motion  of  Atmospheric  Point  Vortices • 

We  consider  the  motion  of  n  atmospheric  point  vortices 
in  the  entire  x-y  plane  without  solid  boundaries.   First, 
consider  a  single  vortex  at  the  origin  (see  Fig.  1).   The 
stream  function  ^     of  this  vortex  of  strength  y  is  the 
solution  (which  vanishes  at  infinity)  of 

(A  -  *<^)^  =  Y  •  5(x)  •  5(y)  (I6a) 

where      5(x)    •    5(y)    =  5(x,y)      is    the   txiro-dimensional  Dirac 
6-function;    thus, 

vf  =  U  K^(Kr)  (16b) 

where   K  (Kr)   is  the  modified  Bessel  function  of  the  second 


2 


*  y2)l/2. 


kind  and  r  =  (x   +  y  )    •   The  strength  y  is  obtained  by 
integrating  equation  (l6a)  over  an  arbitrary  region  /^  in- 
cluding the  vortex;  thus, 

Y  •  j/  5(x)  .  6(y)dx  dy  =  [J  {/\   -  K^)^   dx  dy       (l6c) 

n  k 


or. 


r  =    Y    vf  ds  -  K^  J  J  \lf  dx  dy  (l6d) 

where  the  subscript   v  denotes  differentiation  in  the  direc- 
tion of  the  outward-drawn  normal  to  the  contour  C,   The  line 
integral  is  the  circulation;  and  the  sign  of   y  has  been 
chosen  such  that  positive   y  corresponds  to  a  counter  clockwise- 
rotating  vortex  (cyclonic).   Since  i/\  -  K    )i|/  =  0  everywhere 
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except  at  the  origin,  the  vortlcity  distribution  of  an  atmos- 
pheric vortex  is  given  by 

2 
^(1)  =  ^^  =  K^^l/  =  -  1^  K^(Xr).  (I6e) 

We  note  that  the  vorticity  ^^    '      and  stream  function  ^     have 
the  same  sign.   The  velocity  (tangential)  distribution  is 

^r  =i^^'o(^^^  =  i^^i^^^^-  ^^^^^ 

Near  the  origin   (r  — >  0) 

which  behaves  like  the  ordinary  logarithmic  point  vortex. 

At  a  large  distance  from  the  origin   (r  — >  oo  ) 

_1 
vl/  ^  (Kr)  ^  •  e'^^  (I6h) 

which  damps  out  faster  than  the  ordinary  point  vortex.   We 
see  that  the  atmospheric  point  vortex  defined  by  equations 
(16)  has  the  qualitative  features  of  a  closed  high  or  loi-r 
pressure  system  in  the  atmosphere,  for  example,  a  hurricane 
(or  typhoon).   The  approximate  representation  of  a  hurricane  by 
an  atmospheric  vortex  breaks  dovm  completely  in  the  immediate 
vicinity  of  the  vortex  point,  corresponding  to  the  eye  of  the 
hurricane;  but  this  deviation  is  to  be  expected  even  ^^rith  more 
realistic  solutions  of  the  Geostrophic  Conservation  Equation 
(11a)  since  both  the  hydrostatic  and  geostrophic  approximations 
become  invalid  there.   The  atmospheric  point  vortex  is  the 
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simplest-possible  closed  model  of  geostrophlc  flow  and  differs 
from  the  ordinary  point  vortex  in  several  respects:   1)  the 
direction  of  rotation  is  coupled  to  the  sign  of  the  perturbed 
depth  h^"*"'   which  corresponds  to  the  perturbed  pressure, 
2)  the  motion  is  rotational  as  expressed  by  equation  (l6e), 
and  3)  comparing  vortices  of  the  same  strength,  the  Influence 
of  an  atmospheric  vortex,  as  expressed  by  equation  (l6b),  is 
less  than  that  of  an  ordinary  vortex. 

The  motion  of  ordinary  point  vortices  have  been  studied 
extensively  since  the  early  work  of  Helraholtz  and  Kirchhoff 
(e.g.,  see  Lin  1914-3  and  Lamb),  However,  all  the  analyses  with 
which  we  are  familiar  rest  on  the  implications  of  Irrotational 
motion  and  Bernoulli's  law  without  direct  reference  to  the 
vorticity  conservation  equation  (II4.C).   In  our  study  of  the 
motion  of  atmospheric  point  vortices,  we  must,  of  necessity, 
start  with  the  Geostrophic  Conservation  Equation  (11a).   This 
derivation,  which,  of  course,  also  holds  for  ordinary  point 
vortices,  is  carried  out  in  Appendix  II-B  by  initially  rcn- 
slderlng  lumped  finite-area  symnietrical  distributions  which 
do  not  overlap  for  all  tlm.e   t  >  0,   In  the  following  we 
give  the  results. 

For  n  atmospheric  point  vortices  at   (x.,y.),  the  stream 

function  ^     satisfies 

n 
(A.  -  ^-^H  =  y~  Yi  •  5(x  -  x^(t;))  .  5(y  -  y^(T))    (17a) 
1=1 

and 
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vKx,y,T)  =  ^,J__   YiK^(  r. )  (17b) 

1=1 

?  2  1/2 

where  r.  =  [ (x  -  x. )   +  (y  -  y- )  ]    •   The  velocity  of  the 

k   vortex  is  obtained  by  differentiating 

n 

1=1 

that  is,  t?he  regular  part  of  equation  (17b)  at  the  vortex 

point   (^i^jy^)*'  ^"^ 

(1)     3   ,    ,      ^s 
^(k)  =  -  di^  ^(k)(^k'yk'^^ 

T  Y^    ^Yi  -  "V".  )  '   I         dx, 
=  I?  Z_  ^1  -^-  •  ■'oC'-ik)  =  dT^  (18b) 

1=1  ^^ 


and 


lA 


^(k)  =1^  ^(k)(^k»yk''^^  ^19°^ 


2n  [_   Yi  T:, •  ^^^o^^'^ik^  =  dT 

1=1        ^'^ 


-1  \   ..   '  k   ^i'  .  ,,«  ,..„   ^  _  ^'Mc 

1=1 
iA 


2  2  1/2 

where   r^^^  =  [  (x^  ~  ^1  ^   ■*"  ^^k  "  ^i  ^  ■'    *   Kirchhoff  has  shown 

that  we  can  express  equations  (l3a,b,c)  in  a  more  elegant  way, 

and  the  motion  of  the  k   vortex  can  be  described  by 

dx,    -s  , 
-^k  dT^  =  %  (19a) 

k 
where   the   Kirclihoff  function 


!-.!^ 


{ -rv:. 
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Equations  (19a, b,c)  are  a  system  of   2n  first  order  ordinary 
non-linear  differential  equations  which  are  more  susceptible 
to  analytical  or  computational  treatment  than  the  Geostrophic 
Conservation  Equation  (11a),  especially  if  the  niimber  of 
vortices  is  small. 

A  related  solution  is  obtained  if  ;^e  transform  the 
Geostrophic  Conservation  Equation  (11a)  to  polar  coordinates 
(r,©,T)   and  seek  solutions  satisfying  the  condition  of  zero 
radial  flow.   Then  equation  (11a)  becomes 

[^(ril/p)^  -  ^^^'\   =  0  (20a) 

with  ^     =   0,3  Thus  we  obtain 

where   P(r)   is  an  arbitrary  function  depending  on  the  initial 
conditions.   If  we  pick  P(r)  =  5(x)  •  5(y),  we  obtain  the 
atmospheric  point  vortex  (unit  strength),  ^   =  -{l/2rr)K    (Kr  ); 
for  an  arbitrary  F(r)   the  solution  of  equation  (20b)  is  given 
by  equation  (II-5b)  in  Appendix  II-B.   If  P(r)   is  a  smooth 
symmetrical  function,  then  so  is  the  corresponding  solution   vj/. 

The  motion  of  a  single  atmospheric  vortex  Imbedded  in  a 
continuous  field  is  of  some  interest,  particularly  in  view 
of  possible  application  to  the  problem  of  predicting  the  motion 
of  a  hurricane ^s^  We  write  the  stream  function  ^     as  the  sum 


3k.* 


of  an  axlally-syinmetric   fionction     ^   ,   which  may  be   sihgular> 


and  a  continuous   field     \|f. 


^  =  ^o^^o^  "^  \l'l(x,y,'T) 


2  _  r..    ..  /_^^2  ,  r__    ._  ,^^^2 


(21a) 


where  r'^  =  [x  -  x^i'v)]''  +  [j  -  j^{t:)]^     and  vjf^  is  the 
pat-ticular  solution  of 


(A  -  K    )^   =  F  (^o) 


O'  O' 


(21ta) 


Substituting  equation  (21a)  in  the  geostrophic  conservation 
equation  (11)  we  obtain 


o  L  ■■' 

(22a) 


where 


u^^^  =  -{^^  +   il/-,  )y   and   v^^^  =  (\1/q  +  vl/^)^* 
Following  the  same  argument  as  in  Appendix  II-B,  -vje   take 


(22b) 


and 


Fo(r^)  =  5(r^)  =  6(x-x^)  •  0(7-7^) 


dx. 


dT     (o)     '^1  y 


5^-  v(l^ 

dT   ~  ^(o) 


(^'iV. 


(23a) 
(23b) 


where  {u}    {,    v,~{)      is  the  velocity  evaluated  at  the  vortex 
point   (x  ,y  )   using  the  regular  part  of  \|/,  namely  \1/,  ,   Then 
the  expression  within  the  bracket  in  equation  (22a)  vanishes 
and  vre  obtain 


.  f.-?' 
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d^^  (A  -  '<^)H'i  =  0.  (21+) 

Equations  (23a, b)  and  (2i+)  are  a  system  of  three  equations 
to  be  solved  for  ^1/-,  (x,y,T) ,  x^(t)   and  j^{i),    given  the 
initial  conditions.   Although  we  must  take  ^     =  ■•{Y/2'n)K    (kt    ) 
in  this  formulation,  we  believe  that  a  more  realistic  choice 
of  ^        (such  as  a  continuous  axially-symraetric  fxmction)  may 
be  used  to  good  approximation,  especially  if  ^       is  a  rather 
sharply -peaked  function  as  compared  to   y];,  , 

The  direct  applicability  of  our  concept  of  atmospheric 
vortices  to  meteorological  motions  depends  on  the,  as  yet, 
tentative  assumption  that  we  can  interpret  the  depth  h   of 
the  atmosphere  at  rest  such  that  the  perturbed  depth  h 
(or  \|/)   satisfies   (^  -  K  )\1'  =  0  approximately  over  large 
areas  of  the  region  of  interest;  this  condition  may  be  overly 
restrictive,   V/e  also  point  out  that  flov/s  with  large  curvature 
such  as  near  the  eye  of  a  hurricane  deviate  considerably  from 
geostrophic  flow  so  that  we  certainly  cannot  expect  good  quan- 
titative agreement  there.   However,  we  believe  the  general 
techniques  and  concepts  presented  in  this  paper  will  be  help- 
ful to  meteorologists* 
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Appendix  II-A,  A  Linear  Perturbation  on  a  VJest-East  Flovj  by 
the  Ordering  Approximation  Procedure, 

A  linearization,  related  to  the  Geostrophic  Conservation 
Equation  (11)  and  its  linearized  form  equation  (13  )>  can  be 
obtained  by  the  ordering  approximation  procedure  in  a  manner 
similar  to  the  derivation  in  Appendix  I-B,  but  relaxing  all 
conditions  1),  2)  and  3)  somewhat.   Consider  the  resulting 
stretching  transformation  with  respect  to  a  small  parameter  a, 

^  =  ax,      t;  =  at  (II-l) 

with  the  remaining  variables  unstretched  and  the  associated 
perturbation  expansion  on  a  uniform  west-east  flo^^r  u  : 

h  =  h^(y)  +  ah^^^?,y,'i:)  +  a^h^^^^,y,T)  +  ... 

u  =  U  +  au^^^(4,y,T)  +  a^u^2^^,y,T)  +  ...   (II-2) 

V  =  av^^^^,y,T:)  +  a^v^^^  (?,7,t)  +  ...   . 

Putting  (II-l)  and  (II-2)  into  equation  (lb)  and  (2a,b),  the 
solution   (h^^^,u^^^,v^-^)   satisfies 

With  the  geostrophic  conditions  on  the  velocity  coijponent  u 
only 

U  =  -g(h^)  /f,   u^l^  =  -ghl^Vf    (II-3b) 
and  v^  '  =  0,  which  here  does  not   imply  h^   =  0;  f  =  constant. 
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Appendix  II-B,   Derivation  of  the  Equations  of  Motion  of 

Atmospheric  Point  Vortices oNit' 

Consider  each  of  n  vortex  points   P.   covered  by  a  smooth 

S3mmietrical  distribution  function  \^  ^(^^^  )j  ^w  ~ 

1      1 

[(x-x.(t:))^  +  (y-y.(T))  ]  '  ,  over  a  finite  circular  area 

a.,    with  F  =  0   outside   /^^   (Fig.  II-l).   It  is  implied 

that  these   P.   lumped  distributions  do  not  overlap  for  all 

t;  >  0,   Then  the  solution  of 

n 

(A  -  >^^)^  =  Yl   ^i^^'^xx.  ^  ^"-^^ 

1=1        ^ 


n 

2=1 


n 
=  II^i5(^xx,^  (=^^-^^^ 


where 


is  also  a  smooth  symmetrical  function  with  respect  to  each  P.j 
y.      is  the  vortex  strength  (constant),   3y  the  Geos trophic 
Conservation  Equation  (11a)  and  equation  (II-I4.)  we  v/rite  form- 
ally 

'^  i=l 


d^illYiFCr^,.)   =0  (II-6a) 


or 
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n  Y..F  (r^^  ) 


(II-6b) 
dx.         dy.  ] 

Since  by  construction  F  (r  ,  )  =  0  for  all  points  outside  of 
K?    equation  (II-6b)  is  identically  satisfied  there;  and  we 
need  only  consider   (x,y)   inside   ^j^,  say,  for  the  k  "  vortex 
point   P,  .  Then  equation  (II-6b)  becomes 


([(x-xj^)u(l^  ^  (y-yk)v(l^] 


""^^k    ■  (II-7) 

^nr  ^^k  ] 

-  t(--k)  d^-^  (^-^k)  dFrV  =°- 

Now,  by  equations  (llb,c)  and  (II-5a)  we  can  write 

n 

i=l         ^ 

But  since   $(r   )   is  a  symmetrical  function  about  the  point 
xx^ 

P,.  independent  of  Q,      we  have 
{(y-y^)|e  -  (x-x^)|y}  i(r^J  =  -  ^   S(r^^)  =  0.     (IT-9) 

Hence  in  equation  (II-7)  we  can  replace   (u^   ,v''   )   by 
^^(k)»^(kr   where 

•   uj?-!  =  -  I-  ^n  ^,      v!?"!  =  —  \\if,   X  (IT-lOa) 

( k )     oy  ^  ( k ) »    ( k )    dx  M  k ) 


and  ^(y)      is  the  "regular  part"  of   \|/: 
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♦(10  'i-^i^'^^xx.'-  '"-lO"' 

1=1  ■'■ 

lA 


Thus  equation  (II-7)  becomes 


^k^'(^xx,  ^O-(k)  -  H^^^°^  «k  -^  (^(k)  -  a^^^i-  ^k}  =  0-  (^^-11^ 


If  we  now  choose     F     to  be    such  a  fiinctlon  that  11m     F(r        ) 


R  _^o   ^^k 


=  5(x-x,  ,y-y,  )  =  6{x-x,  )•  5(y-y,  ),  the  Dlrac  5-f unction,  then 
equation  (11-11)  Is  satisfied  by 

u!?-!  -  XT^  ==  0  (II-12a) 

(k)    dT 

v|Jj-^=0  (II.12b) 

where   (^n^^^v)-  {)   Is  evaluated  at   (x,y)  =  (x  ,y  ),   The 

symmetry  of  the  5-functlon  Is  Implied  by  the  symmetry  of 

f(r  ^  )   which  by  equation  (Il-ii.)  and  (II-5)  becomes 
^k 

$(r  „  )  =  -(1/2it)K  (K:r    )»   Thus  by  choosing  the  distri- 
k  ^^^k 

bution  function  P  properly  we  have  replaced  the  Geostrophic 
Conservation  Equation  (11a)  by  equations  (II-12a,b)  describing 
the  motion  of  atmospheric  point  vortices. 


FOOTNOTES 

\1/   This  paper  is  based  on  work  cosponsored  by  the  Office  of 
Naval  Research  and  the  Geophysics  Research  Directorate, 
Air  Force  Cambridge  Research  Center,  under  contract 
N6orl-201,  T,  0.  No.  1  with  the  Office  of  Naval  Research, 
Reproduction  in  vjhole  or  in  part  permitted  for  any  purpose 
of  the  Ue  S,  Government, 

\2j'      Except  near  the  equator  ^^fhere  the  horizontal  component  of 
the  Coriolis  force  approaches  zero, 

\3/   Just  such  a  derivation  has  been  attempted  in  an  unpub- 
lished paper  by  J,  B,  Keller  and  Lu  Ting  at  N.Y.U,  with 
some  success, 

\k/  Stoker  (195>3)  has  considered  such  systems  in  his  study 
of  frontal  actions  in  meteorology, 

\^     Actually,  this  condition  is   v^   f   =0;  but  we  consider 
v     not  identically  zero. 

\6/  Work  on  higher  order  approximations  have  been  reported 
recently  by  Charney  (1955),  Bolln  (1955)  and  Thompson 
(unpublished  manuscript). 

\^  With  A,  Troesch,  Institute  of  Mathematical  Sciences,  N.Y.U, 

\8/  Like  Prandtl's  boundary  layer  equations, 

\^  For  economy  of  notation  Xire  make  this  deviation  from  the 
main  body  of  the  paper, 

\10/  This  result  is  in  disagreement  with  conclusions  reached 
by  Bolln  (1955)  for  a  tvjo-layer  model  using  a  less  con- 
sistent procedure. 
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sll/  Suggested  by  E,  Isaacson, 

sl2^  Charney  informs  us  that  the  more  recent  numerical  pre- 
dictions have  been  made  using  a  comparable  computational 
procedure  (private  communication). 

\13/  A  series  of  interesting  papers  on  the  apparently  successful 
prediction  of  hurricane  motion  by  Japanese  meteorologists 
has  been  brought  to  our  attention,  e.g.,  Sasaki  and 
Miyakoda  (1955 )•   ^'^e  seem  to  have  some  common  general 
viewpoints  but  their  approach  and  method  are  on  a  more 
physical  basis. 

\lli^  With  Bo  Z^amino,  Institute  of  Mathematical  Sciences  and 
Department  of  Physics  (University  Heights),  N.Y.U. 

nJ.5/  We  use  a  more  descriptive  notation  here  for  the  distance 
r„„  E  r.   (cf.  Section  i;). 
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